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STOCHASTIC EXCITATION OF STELLAR OSCILLATIONS 

R. Samadi 1 

Abstract 

Excitation of solar oscillations is attribued to turbulent motions in the solar convective zone. 
It is also currently believed that oscillations of low massive stars (M < 2 Mq) - which possess 
an upper convective zone - are stochastically excited by turbulent convection in their outer 
layers. 



A recent theoretical work (Samadi and Goupil, 2001; Samadi et al., 2001) supplements and 



reinforces this theory. This allows the use of any available model of turbulence and emphasizes 
some recent unsolved problems which are brought up by these new theoretical developments. 



1 Introduction 

The excitation of the solar oscillations results from the action of the Reynolds tensor and of the 
turbulent fluctuations of entropy; the latter comes from a heat exchange between the oscillations 
and the turbulent elements. The nature and the properties of the source term of excitation due to 



the Reynolds tensor are well established (Goldreich and Keeley, 1977), although the evaluation of 



this term remained crudely approximated. On the other hand, the source of excitation which takes 
its origin in the fluctuations of the entropy, had given rise to major controversies (e.g. poldreich 



et al., 1994, GMK hereafter). 



In Samadi and Goupil (2001, Paper I hereafter), inconsistencies in the available theories were 
removed which led to a formulation which generalizes results from previous works, and is built so 
as to enable consistent investigations of various possible spatial and temporal spectra of stellar 
turbulent convection. It also shows that the actual entropy source term results from the advection 
of the turbulent fluctuations of entropy by the turbulent motions. These results are summarized 
in Section and Section ||. 

The unavoidable free parameters introduced in the formulation are adjusted in order to obtain 
the best fit to the solar seismic observations. In a second step, the comparison between the 
computed amplitudes and the observations as well as the use of theoretical arguments allow us to 
determine the ingredients in the theory which are still defective. Section 4 then lists some unsolved 
problems which are brought up by our improvment of the theoretical approach in Paper I. Final 
comments conclude about elements of the theory which still remain to be developed. 

2 Theory of the stochastic excitation 

The stochastic excitation mechanism results from a forcing of the stellar material by the turbu- 
lent motions and by a heat exchange due to the turbulent fluctuations of the entropy; in other 
words, the acoustic power generated by the turbulence excites resonant modes of the stellar cavity 
(oscillations). 
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In the presence of turbulence one shows that the radial oscillation velocity v osc obeys the 
inhomogeneous wave equation (see Paper I for details): 



d 2 A,- d 

Po I 

with 



P0 1 dt 2 



L^j [v osc ] + V{v osc ) = ^(ft + V/i t ) (2.1) 



ft = -V : ( Po uu) (2.2) 

d - - d 

— Vh t = -V ( a s —Ss t - V.(a s us t ) + stu.Vcts ) (2.3) 

the turbulent Reynolds stress and the source term due to the turbulent entropy fluctuations 



respectively. The expression for the wave operator L is given in Paper I. In Eqs^lO , 2/2, 2.3), 
p, s, c s , g and u respectively denote density, entropy, sound speed, gravitational acceleration and 
velocity of the turbulent elements; the subscript refers to average mean quantities. 

The first term in the RHS of Eq.( [2.3| ) is due to the Lagrangian entropy fluctuation 6s t- The 
last two terms are due to the buoyancy force associated with the Eulerian entropy fluctuations 
St- These terms are found to contribute to the excitation as much as the Reynolds source term. 



Other source terms in the RHS of Eq.(2T) are found negligible in Paper I and are discarded here. 

The operator T> in Eq.(2.1) is responsible for generating both a dynamical damping and a 
shift of the oscillation frequency due to the action of the turbulent elements. Contribution of 
the operator T> to the damping is expected to be small compared to those of the other damping 
processes, and is assumed to be included in the global damping rate 77 which takes into account 
all damping processes. 

Assuming no turbulence (u = 0) the velocity field is simply related to the pulsational radial 
displacement as v — ~^£,r [?■> t) where <^(r, t) is the undriven pulsational radial displacement in the 
absence of turbulence. Eigensolutions of Eq.( |2.l| ) complemented with boundary conditions, can 
be written, in the absence of turbulence, in terms of the displacement £ r (r, t) = £ r (r) e"*" * where 
u)q is the oscillation frequency and £ r (r) is the adiabatic (real) displacement eigenvector which 
satisfies L(£ r (r)) = — a;Q^ r (r). 

When turbulence is present, the pulsational displacement and velocity are written in terms of 
the above adiabatic solution £ r (r, t) and an instantaneous amplitude A(t). Accordingly 

8r osc = l - (A(t) e^ ot + cc) and v osc = ~(-iLJ A{t) £(r) e' wat + cc) (2.4) 

where cc means complex conjugate. The mode energy, averaged over a time scale smaller than 
the mode lifetime i]^ 1 and larger than the eddies time scale, is given by 

,M 1 

E = J a dm (v 2 osc ) = - <L4| 2 > W (2.5) 

where 1=1 Pa d 3 x (£* .£) r is the mode inertia and (|^4| 2 ) is the mean square amplitude. Finally 
Jv 

a general expression for the power, P, going into each mode is established in Paper I: 

p.f-,nv4 w„h < M n-j-^cs+ts) ( 2 .6) 

where C\ and C| are the turbulent Reynolds stress and the entropy contributions respectively 
given by: 

C\ = d 3 x dre-^ \ d 3 r((p u j u i Vi£) 1 (p u j u i V j £) 2 ) (2.7) 



— OO 



C 2 S = Id 3 



3 x dTe-**> T / d 3 r^/i t V.er) x (h t V.ir) 2 ) (2.8) 
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where xq is the average position where the power is evaluated whereas f and r are related to the 
local turbulence. Subscripts 1 and 2 refer to the spatial and temporal positions [xq — %, — ■k] and 
[x + L J] respectively. 

2.1 Reynolds stress contribution 



It is shown that the terms (pqV^^) and (poV^Cr) m Eq.(2.7) do not change over the length scale 



of the eddies. Consequently, in Eq.(2.7), integrations over r and r only involve the phase term 
e -iui r an( j fourth-order velocity correlations (u^u'jU'i'u'^} where v! = u(xo — r/2, — r/2) and 
u" = u(x Q +r/2,+r/2). 



The Quasi-Normal Approximation (Lesieur, 1997, Chap VII-2, QNA hereafter) is adopted so 



that the fourth-order velocity correlations can be written as a product of second-order velocity 
correlations (ii^u"). The properties of {u^ii") are well known in the Fourier domain (k,Lu) where 
k and uj are the wavenumber and the frequency associated with a turbulent element. The second- 
order velocity correlations can be expr essed formally in term of E(k,uj) the turbulent kinetic 



energy spectrum. Following Btcin (1967 ), the velocity energy spectrum E(k,uj) is written as 



E(k,u) = E(k)xk(u) (2.9) 

where Xk(^) is a frequency-dependent factor which can be related to the time correlation of the 
eddies in the frequency space. Finally a general expression for involving E(k) and Xfc(^) is 
given in Paper I: 

C\ = 47r 3 S J\mp (^\ ' dk^-Xk(u>o) (2.10) 
where Q is a constant anisotropic factor. 



2.2 Contribution of entropy fluctuations 



We assume in Paper I that dst and s t in Eq.(2.3) act as passive scalars. Let E s (k,to) be the 



spectrum of the scalar correlation product (s' t s"). It is also assumed that the frequency-dependent 
component of (sjs") is the same than those of the velocity field correlation product (u^m") and 
that E s (k,Lu) can be decomposed as E(k,uj) in Eq.(|2.9|). 

As a consequence of the above assumptions, it is then demonstrated that the correlation prod- 



uct of the Lagrangian entropy fluctuations {Ss' t Ss") vanishes after integration over r in Eq.^ 
This result may be explained as follows: turbulence and oscillation are coupled through the phase 
term e~ lUaT and through the turbulent time spectrum Xk(u>). Therefore this coupling occurs at 
frequencies close to the oscillation frequency ujq. It then follows that the coupling between tur- 
bulence and oscillation involves eddies of wavenumber k 3> k osc where k osc is the wavenumber of 
the mode and k that of an eddy. On the other hand the spatial component of {5s' t 5s' t ') in the 
Fourier space favors eddies with the largest size (k — > 0). These two opposite effects clearly are 
incompatible and lead to a vanishing contribution for the Lagrangian entropy fluctuation. 

This does not happen for the contribution of the Reynolds source term which involves the 
fourth-order velocity correlation product. According to the QNA this term can be decomposed 
in term of a product of two second-order velocity correlations. Coupling with the oscillation then 
becomes non-linear and leads to an effective non zero contribution. Thus only non-linear terms, 
with respect to the fluctuations, can contribute to mode excitation while linear terms do not. This 
may be considered as a general result. 

With an analogous procedure than for the Reynolds stress contribution, we show that the 
entropy contribution can be expressed as: 

/' ('.. <l^ \\. f^E s (k)E(k) ' - 



C s = —jT~ J d x ° [ as ~dr J 9r j dk ~k~2 J dtvxkiLOQ + Lu)xk(u) (2.11) 

with g r (lj,r, f) a function which expression is given in Paper I and TL a constant anisotropic factor 
similar to Q. 
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3 Results in the Solar case 



3. 1 Models for the solar turbulence 

On one hand, the turbulence theory tells us that E(k) follows the Kolmogorov spectrum as E(k) cx 
k p with the slope p = —5/3. On the other hand, observations of the solar granulation allow one 
to determine the turbulent kinetic spectrum E(k) of the Sun. 

Several kinetic turbulent spectra have been suggested by different observations of the solar 
granulation: the "Raised Kolmogorov Spectrum " (RKS hereafter) , the "Nesis Kolmogorov Spec- 
trum" (NKS hereafter) and the "Broad Kolmogorov Spectrum" (BKS hereafter). These spectra 
(Figure [l]) obey the Kolmogorov law for k > k where k is the wavenumber at which the turbu- 
lent cascade begins. They differ in the injection region (fc < k ) where they are characterized by 
different values for the slope p. 

The wavenumber ka is unknown but can be related to the mixing length A as fco — 2tt / (/3A) 
where j3 is a free parameter introduced for the arbitrariness of such definition (Paper I) . Moreover 
the definition of the eddy correlation time , which enters the description of the turbulent excitation, 
is somewhat arbitrary and is therefore gauged by introducing an additional free parameter A. 



3.2 Oscillation power 

The power injected into the solar oscillations is related to the rms value v s of the surface velocity as 
v^(ujq) — £, 2 (r s ) -P(wo) / 2r]I where r s is the radius at which oscillations are measured. Observations 
of the solar oscillations provide v s and 77 such that P(u>o) can be evalua ted f rom pbservatio ns. On 



the theoretical side, the power P(u>q) is computed according to Eqs(2J3, 2.6, 2.1C , 2.11) for a 
calibrated solar model. The physical ingredients of the model are detailed in Bamadi et al. (2001 , 
Paper II hereafter). As was done in Paper II, the free parameters ((3, A ) are adjusted in order to 
obtain the best fit to the solar seismic observations by [Libbrccht (1988| ) 

Theoretical arguments support the fact that the current theory of stochastic excitation is less 



justified for the low frequency modes than for the high frequency ones (see below section 4.1). 
Therefore in the following (3 and A are adjusted in order to reproduce, as best as possible, the 
frequency dependency of v s (u>o) at high frequencies (y > 3.5 mHz). 

Results of the fitting are shown in Fig. [|. As a consequence of the adjustement, all the kinetic 
turbulent spectra fit well the solar observations at high frequency [y > 3.5 mHz) while the main 
differences are observed at low frequency. The overall best agreement is obtained with the NKS. 

Moreover we find that the entropy source term significantly contributes to the excitation pro- 
cess of solar-like oscillations by turbulent convection, in agreement with the results of GMK and 



simulations by Stein and Nordlund (1991) 



4 Some open issues 



4.1 Static and dynamic properties of the turbulent medium 



As stressed by Rieutord et al. (200C ), differences at low wavenumber between the kinetic spec- 
tra obtained from different observations of the solar granulation are a consequence of uncon- 
trolled data-averaging procedures. Mean average (static) properties of the kinetic spectrum at low 
wavenumber (mesogranul ation) are thus not yet well represented by the observations of the solar 
surface. In the same way, Nordlund et al. (1997 ) demonstrated with a numerical simulation that 
observations of the solar granulation cannot provide a meaningful determination of the turbulent 
properties. Therefore, the mean average (static) properties of the different turbulent spectra - 
suggested by the current observations of the solar gran ulation - a re no t yet well established. 

On t he theoretical side, the derivation of Eq. ( |2.10 ) and Eq.( [2.11 ) uses the decomposition of 
Eq.(2.9). E(k) depicts the static properties of the turbulent medium whereas E(k, u>) characterizes 
the dynamic properties. In Paper II different forms for Xk( w ) _ which were suggested by Musielak 
et al. (1994| ) - were tested. From the comparison with the solar seismic data we cannot clearly 



discriminate between the different forms and we finally adopted the gaussian form (Paper II). 
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Fig. 1. Left panel: Kinetic turbulent spectra versus wavenumber k. Right panel: Computed surface 
velocity v s assuming the turbulent spectra plotted on the left panel. A and ft values result from fitting the 
computed v a to the solar seismic observations by Libbrecht (1985) using the observed damping rate r\. 



However it is suggested in [Bamadi (200CD that the decomposition of Eq.( |2.9D is not valid at 
small wavenumber (in the injection region, i.e. for k < feo). Indeed the high frequency modes 
involve small size turbulent elements (k > fc ). At this scale (the inertial range) the lifetime 
Tfc of the turbulent elements are shorter than the timescale ta at which energy is injected into 
the turbulent cascade. Therefore at small scales injection of energy into the cascade appears 
stationary and the dynamic properties of E(k,u/) are only fixed by the lifetime of the eddies. 
The line-width of the function Xk(u) provides a statistical measure of Tfc. As a consequence, the 
decomposition of Eq.(2.9) should be valid for k > k . 

In contrast at large scales, ( k < fco), the lifetime of the turbulent elements is of the same order 
than ta- Therefore the energy transfer between large and small elements is no longer stationary 
and consequently the decompostion of Eq. ( |2.SD is probably no longer valid. 

This statement explains the large changes observed at low frequency between the different 
spectra (Fig. |l|). Indeed the low frequency modes involve the largest eddies and the changes at 
low frequency are therefore due to the differences in the static properties of the turbulent spectra 
in the injection region (see Fig. |). The best agreement is obtained with the NKS which however 
appears as the less realistic spectrum. In contrast the more plausible spectra (RKS and BKS) 
significantly overestimate the amplitudes at low frequency. These results prove that the dynamic 
properties of E(k, u>) is not correctly modeled through the decomposition of Eq.(|2.9[) in the range 
k < kg. Moreover, the decomposition of Eq. ( p^9| ) is also used for the turbulent spectrum E s {k,u>) 
of the entropy fluctuations. Therefore the above conclusions also apply to E s (k,to). 



4.2 Excitation of the non-radial oscillations 



The present theory of stochastic excitation only concerns the radial p-modes. It is shown in 3amadi 
|(2000| ) that the stratification does not affect the transfer of acoutic power into the radial modes. 
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The high £ degree modes propagate in both radial and horizontal directions. The acoustic 
power injected into high £ degree modes is therefore more sensitive to the stratification. 

The current theory is mainly based on the approximation that the stratification and the sources 



of excitation are well decoupled. Stein (1967 ) showed that the acoustic emission arising from the 
stratification is negligible compared to the Reynolds stress emission. However some high £ degree 
modes propagate mainly in the horizontal direction in the excitation region. For such modes it is 
important to evaluate the contribution of the stratification to the acoutic power emission. 

5 Conclusion 



The present paper summarizes the theory of solar oscillation excitation at it is exposed in 3amadi 



and Goupil (2001). In particular assumptions and approximations which were used are recalled 
here with some details. Open issues in this field are also discussed such as the way the turbulence 
is modeled and the case of non-radial modes: 

The way the turbulence is currently modeled at large scale length introduces large uncertainties 
in the computation of the power injected into the low frequency modes {v < 3.3 mHz in the solar 
case). A more robust modeling of the large scales is thus needed in order to extend to low 
frequencies the validity of the current theory. This can be achieved with the use of 3D simulations 
of the solar convection zone (work in progress). 

Current theories only concern the excitation of radial modes. For non-radial modes effect of 
the stratification upon the excitation may be important whereas it is negligible for radial ones. 
Therefore effects of the stratification to the power emission in the high £ degree modes, as well 
as the £ dependency of the current theory of stochastic excitation, should be addressed (work in 
progress). 
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